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Abstract 

We explore the elementary observation that a Markov chain with values in a hnite space M 
with \M\ = m, m > 2, has many different extensions to a compatible n-point Markov chain in 
M”, for all 1 < n < TO. Embedding this phenomenon into the context of stochastic Levy flows 
of diffeomorphisms in Euclidean spaces, we introduce the notion of an n-point bifurcation of a 
stochastic flow. Roughly speaking a n-point bifurcation takes place, when a small perturbation 
of the stochastic flow does not change the characteristics at lower level fc-point motions, k < n, 
but does change at the level of n-point motion. We illustrate this phenomenon with an example 
of an n-point bifurcation, with n > 3. In addition, we present an algorithm for the detection of 
the precise level of an n-point bifurcation and a combinatorial formula for the dimension of the 
vector space of compatible extensions for flows of mappings on M. 

Keywords: stochastic dynamics, stochastic bifurcation, n-point process, stochastic Levy flow, 
Markovian process, algorithmic detection of bifurcations, Stratonovich SDE, canonical Marcus 
SDE. 

MSC2010 subject classification: 37H20, 60J10 (93E03). 

1 Introduction 

It is well-known for Brownian flows that the law of the stochastic flow is uniquely determined by 
the laws of the family of the corresponding 1-point and 2-point motions, see Kunita m- This 
behavior is due to the Gaussian nature of the marginal laws and can be read from the structure 
of the infinitesimal generator. General Levy flows do not exhibit this behavior. In this article we 
construct a simple example, which motivates us to coin the notion of a general stochastic re-point 
bifurcation. 

We explore the elementary observation in the context of stochastic flows that a Markov chain 
with values in a finite space M = {!,... ,rn}, with m > 1, has many different extensions to a 
compatible re-point Markov chain in M”, for all 1 < re < rre. Each extended process in is called 
an re-point lift of the original 1-point motion in M. Gompatibility here means that any projection 
of the process in into an embedded M^, with k < n leads to a Markov chain in with the 
same law, such that, eventually, after any (re — 1) projections, the original law of the 1-point motion 
in M is recovered. Two different lifts to re-points can still have laws that coincide in the A:-point 
process with 1 < A; < re. This means that the laws of these two lifts can only be distinguished one 
from another if an observer studies the statistics of the j-point process for j > k sufficiently large. 
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We recall that given a probability space (0, J”, P), a (time-homogeneous, measurable) stochastic 
flow in a Polish space is a family of parametrized measurable mappings ip = : A4 —>■ 

Ad I t ^ 0, a; € n} such that <po(^) = Idjvi for all cu G n and (ps+ti^,-) = ips{0tuj,-) o ip^^co,-) 
for all 0 ^ s ^ t and w G id. The time variable t here is considered either in N U {0} or in R>o 
and the shift : Id —>■ fl is an ergodic measurable preserving transformation. The n-point motion 
induced by a stochastic flow ip in Ad” is given by (xi,... ,x„) i-7> {ipt{xi),... ,ipt{xn))- A simple 
example of stochastic flow with discrete time is the composition of a random i.i.d. sequence (^n)n>o 
of mappings in A4. Then, for a positive integer t, the flow ipt = ^t{^) ° ■ ■ ■ ? 2 (w) o 

In the context of discrete Markovian dynamics in M, there are three standard approaches on 
studying the extensions to an n-point version. The first approach consists in restricting to the 
random dynamics generated by composition of sequences of random i.i.d. permutations in M 
(stochastic flow of bijections). This approach, although quite restrictive, is particularly interesting 
since it allows easily to embed the results from the finite space M into a flow of diffeomorphisms in 
an Euclidean space: just use an appropriate choice of the vector fields and Levy noise involved in a 
Stratonovich SDE for general semimartingales of Markus type, as in Kurtz, Pardoux and Protter 
P (see Section n. The second approach is again a stochastic flow, but here one allows for the 
composition of random mappings from M to M, which permits the coalescence of particles. As 
in the bijective approach, the diagonal is still invariant (in fact, an attractor here). Einally, the 
third approach consists in considering a Markovian lift independently of any flow in M, hence one 
misses many symmetries in the transition probabilities in M^, moreover, (sub-)diagonals may no 
longer be invariant. Think, for instance, of the simplest example: M = {1,2} with all possible 
transition probabilities equal to 1/2 and a lifted process in = {(1,1), (1,2), (2,1), (2,2)} with 
all the transition probabilities given by 1/4. This dynamics in is clearly not generated by a 
flow of mappings. 

In this article we focus on the dynamics generated by stochastic flows, i.e. in the first and 
second approaches. The main motivation of this choice is that it allows us to introduce the notion 
of an n-point bifurcation illustrated by the example of Section 2 and 3. 


1.1 The notion of an n-point bifnrcation 

Given invariant measures for the n-point process, the invariant measures in lower /c-point processes, 
k < n, can be recovered using the following well expected result: 

Proposition 1. For 1 ^ k ^ n consider the projection 

pI : M” (xi,...,Xn) (xi,..., x^,..., x^), 

and let p he an invariant measure of a Markov process X in M”. Ifp^{X) is also a Markov process, 
then the induced measure {p^)*pi is an invariant measure for the process p^{X) in 

Proof. Eor convenience we drop the superscript n whenever possible. Let P{x,A) be the family of 
transition probabilities of the process X in M” for x G M” and subsets A C M”. The fact that 
the projection Pk{X) generates a Markov process in means that the transition probabilities 

in denoted by P(pk{x), B), is well defined for any B C and it is given by 

P{pk{x),B)=P{x,p^\B)) 
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for all X G M"". Now, by the induced measure theorem: 




j Pix,PkHB)) dp{x) 

M" 

J P{Pk\y)^Pk\P)) d{pk)*p{y) 
I P{y,B) d{pk)*p{y). 


□ 


The notion of an n-point bifurcation: Classically, in dynamical systems, a bifurcation occurs 
when, changing a parameter of the flow, the topology of the support of an invariant measure 
is affected somehow, normally splitting into two or more disjoint invariant domains. This is well 
known and well studied for deterministic dynamical systems, where the precise dehnition is based on 
breaking local topological equivalences of the flows (among many classical references see, e.g. Katok 
and Hasselblatt [8] and references therein). For stochastic systems generated by Ito-Stratonovich 
equations, the bifurcation is mostly considered as a change of sign of the Lyapunov exponent, see 
e.g. L. Arnold m, 13]. These two situations have in common the fact that they are observing 
a breaking in the topology of the support of invariant measures, but at different levels: In the 
deterministic case, the invariant measures are considered in M^, with trivial extension to M” as 
the respective product measure; in the stochastic case, the sign of the Lyapunov exponents points 
to properties of the invariant measures in M^. See the explicit example by Baxendale [1], where 
a bifurcation happens for Brownian motions in the torus: i.e. the top Lyapunov exponent change 
the sign but the law of the 1-point motion is not affected, see also [5]. 

Going further from these 1 and 2-point phenomena to re-point motions, we propose the following 
natural extension of the definition of bifurcation for more general stochastic flows in Polish spaces: 

Definition 1. Let be a family of stochastic flows acting on a Polish space A4 indexed by a 
parameter e with values in an interval. 

We say that exhibits an n-point bifurcation, n 2, at a parameter value Eq if it satisfies the 
following 

1. The mapping e p^ is continuous with respect to the weak topology at all values e except Eq. 

2. There exists an invariant measure with respect to p‘^° such that tor any e > Eq there exists 
an invariant measure pp with respect to p^ on Ad"' satisfying that 

supp{p^°) is not homeomorphic to supp{if), 
and there exists a projection pf: defined in Proposition [I] such that 

supp{{pl)iP°) = supp{{pl)if). 


□ 
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A first example is given in Section 2. Since each invariant measure on can have many 
lifts to invariant measures in higher levels these lifts can exhibit more than one bifurcation, 
at different levels k. Moreover, the same bifurcation on A:-point motion can have projections into 
different invariant measures on Al^ (depending on the sequence of projections). This is going to 
be illustrated in Section 3.1. 

As for comparison (following Kunita [TO]), we consider a homogeneous Brownian flow ip in the 
group of diffeomorphisms in an Euclidean space, with the infinitesimal mean 

b{x) = lim \ {E[iph{x)\ — x) , V X G 
h—>-0+ h 

and the infinitesimal covariance 

a(x,y)= lim ^ (E [((^/^(x) - x) ((^^(y) - y)*]) , Vx,yGlR'^ 

/i^0+ fl 

whenever defined. Given certain regularity conditions on this parameters (satisfied for instance by 
flows of SDE generated by smooth vector fields with bounded derivatives) the law of ip in the group 
of diffeomorphisms is determined by a{x,y) and 6(x) |10l Thm. 4.2.5, p. 126]. In other words, the 
law of a (homogeneous) Brownian flow (hence the law of its n-point motion, with n > 2) is fully 
determined by the laws of its 1-point motion and its 2-point motion. This theorem tells us that 
classical stochastic flows for SDE generically do not furnish the richness of flows differing only on 
higher n-point motion, with n > 2, as we are looking for. 

The problem we are addressing here is somehow related to the recent results by dost, Kell 
and Rodrigues [7j, where they study conditions under which the transition probabilities (1-point 
motion) in a manifolds can be represented by families of random maps. In the same article, they 
consider further conditions for regularity and representations by diffeomorphisms. Eor this kind of 
problem, where measurability of the flows is considered we refer to Kifer |9| and Quas m- 

The flows we are interested here are also related to the flow of measurable mappings of Le Jan 
and Raimond m (see also M) in the following sense: their flows are constructed from a family 
of Feller compatible semigroups in C'(M"’), n > 1, which preserves the diagonal. They are also 
constructed based on the observation of the statistics of the n-point motion, for n > 1. 

Flow of hijections in M: The celebrated Birkhoff-von Neumann theorem says that nxn bi-stochastic 
matrices lay in the convex hull generated by the m\ matrices of permutations in {1, 2,... , m}. This 
convex set is called the Birkhoff polytopes Pn- There are several proofs of this theorem in the 
literature. For a simple and elementary proof see e.g. Mirsky m- This theory has many interesting 
application, and although already very studied, still has some good open problems, see e.g. Pak 
HSj. For instance, despite its relevance, for higher dimensional Birkhoff polytopes Pn, there is 
no formulae for the volume of P^, only recently an asymptotic formula was obtained by Canfield 
and McKay |6|. In the context of our article, for random dynamics generated by i.i.d. random 
mappings, it means that a stochastic flow in M is a flow of permutations if and only if the matrices 
of transition probabilities of 1-point motion is bi-stochastic. Moreover, in the Birkhoff polytope 
language, what we are exploring in this article is the fact that, in general, except for elements in the 
wedges of the polytope Pn, the bi-stochastic matrices has a non-unique representation as a linear 
combination of the vertices of Pn (in fact, Pn is contained in a {m — l)^-dimensional subspace and 
has m\ vertices). 

The main objectives of this article are the following: 
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1. To illustrate the n-point bifurcation of Definition [T] beyond the well known examples men¬ 
tioned above (for n = 1 and n = 2). Our example starts with a hnite space M which exhibits 
n-point bifurcation, for n > 3. This bifurcation is then embedded in an Euclidean space as 
a flow of diffeomorphisms generated by a generalized Stratonovich equation (Markus type as 
in m) driven by a Levy noise . 

2. Given two different stochastic flows in a finite state space with the same transition proba¬ 
bilities in 1-point motion, we present an algorithm to establish what is the minimal positive 
integer k, such that the transition probabilities are different, or yet, the minimal positive 
integer k such that the invariant measures have topologically non homeomorphic support. 

3. We present a formula for the dimension of the space of all possible vector space of probabilities 
on the mappings from M to M such that the corresponding stochastic flow (generated by 
composition of i.i.d.) satisfies a prescribed /c-point transition probabilities, with 0 < k < m. 

Items (1),(2) and (3) above are treated in Sections 2, 3 and 4 respectively. 


2 Example of an n-point bifurcation with n > 3 

The purpose of this section is to show an example for the higher order n-point bifurcation, which 
does not change the characteristics (transition probabilities) in lower order A:-points for k < n. 
We construct a discrete stochastic flow in M = {1,2, 3,4, 5,6} generated by the composition of 
i.i.d. random bijections with the following property: When the law of the random bijections are 
perturbed by a parameter e G [0,1] in a certain direction, the invariant measures of the 6-point 
motion have different supports for e = 0 and for e G (0,1] but the laws of the corresponding 1- 
point and 2-point motions remain constant with respect to e. Hence, this family of stochastic flows 
indexed by e exhibits a bifurcation for the n-point motions for 3 ^ n ^ 6. This result is embedded 
into flows in in the next section. 

The basic construction: We adopt the notation /iii 2 i 3 i 4 i 5 i 6 the function 

( 1 , 2 , 3 , 4 , 5 , 6 ) !-)■ 

with ii,... ,im G M. The following double notation is convenient for our example: 

(1,2,3,4,5,6) = ((1,2),(3,4),(5,6)) = (a,6,c) 

and denote by a, b and c the flips of each double entry, for example (a, 6, c) = (2,1,4, 3,5, 6). 
Consider the group 


G — {fabc, fabcj fabci fabci fabci fabci fabci fabc} 


with multiplication given by the composition, and its proper subgroup 


LI — {fabci fabci fabci fabc}' 


The 6-point motion of a flow generated by composition of i.i.d. bijections with law concen¬ 
trated on this subgroup, say 


1 

4 


[^fabc + ^fabc + ^/abc + 


( 2 . 1 ) 
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has random trajectories with the following property: for each initial condition in the cor¬ 
responding random orbit of the process is concentrated on at most 4 points ont of 6® possible 
elements of M®. For an initial condition which does not belong to any snbdiagonal (i.e. such that 
their entries are all different from each other), the support of the invariant measure is concentrated 
on exactly 4 elements. On the other hand, the orbits of elements outside any sub-diagonal of the 
6-point motion generated by the e-perturbation in the law, with e > 0, 

\ [^fabc + + ^fab,] + 2 i^hbc + ^Lb. + ^fabb + ^Lb, “ ^fab. “ ^,1. “ h^b, “ ^Lb,] > (2-2) 

has invariant measures supported on exactly 8 elements. Moreover, one easily checks by inspection 
that, due to appropriate cancellations, the transition probability of jumps from a pair of points to 
any other pair of points does not depend on e, i.e the law in is constant. The same happens 
for the law in M^. 

The splitting on the number of connected components of the support of the invariant measure 
implies that there must exist an n-point bifurcation for 3 < n < 6. In the next section we shall 
construct an algorithm to find out exactly at which level n the bifurcation occurs. 

Embedding the n-point bifurcation into flows of diffeomorphisms: Suppose v \s a. 

probability measure on the set S' = M —>■ M of all mappings of M to itself. Suppose (Cn)nGN are 
i.i.d. in S with distribution v and define 

Rfi ■— Cn ® ^n—1 O • • • O 

That is, Rn is the discrete time random walk of the semigroup S generated by the probability 
measure v. To obtain a continuous time version, consider Tj := where {'Kt)tpo is an independent 
standard Poisson process. For / G S let K[f) denote the orthogonal linear mapping of IR™ to itself, 
which sends the canonical basis vectors ei,... , to ey(i),..., ej(m) respectively. Then Xt = K{Yt) 
is a Levy flow of linear mappings of IR™ to itself. If v is supported on the subset of S consisting 
of invertible mappings then Xt is a Levy flow of diffeomorphisms. For further properties on the 
representation of Markov chains by sequences of random maps we refer to [9] . 

Again, as in the discrete case, the laws of the 1-point motion and 2-point motion does not 
depend on the parameter e. But for e = 0 the support of the invariant probability measure for the 
6 -point motion which includes, say the element (ei, 62 , 63 , 64 , 65 , e^) is given by the set of 4 elements 
(4 connected components) 

{(^/(i)’®/(2)T/(3)>e/(4),e/(5)T/(6))> for all / G Flj, 

which splits, for e > 0 into the 8 connected components 

{(®/(i)’^/(2)T/(3)T/(4),e/(5)!e/(6)), for all / G g|. 

Hence, for this stochastic flow of diffeomorphisms, an n-point bifurcation holds at e = 0 with 
3 < n < 6. This example illustrates the well known fact that the Kunita’s result mentioned in the 
Introduction m Thm. 4.2.3] for stochastic classical Brownian flows does not hold for Levy flows. 

3 Detection of the level of an n-point bifurcation 

As before, consider the finite space M = {1,... ,m} with m ^ 2. The purpose of this section is 
to answer the following question: given two transition probabilities at the m-point motion, whose 
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projections coincide at level of fc-point, k < m. What is the lowest level of projections on n-points, 
k < n < m where they differ? Below, we set up an approprioate algorithm, which we apply to 
show that the example of Section [2] exhibits an n-point bifurcation, for n = 3 and n = 5. 

3.1 The background and the algorithm 

For each 1 ^ n ^ m, let be the m”’ x m” right stochastic matrix (i.e. the sum of rows are 
1), whose entries are the transition probabilities among the elements of M"’, in the lexicographical 
order. Since we are dealing with flows, the transition probabilities of can be obtained from 

the projections p” for any r € {1,... , n} defined in Proposition [T] More precisely, for all 1 ^ r ^ n 
and (ii,..., in), (ji, ■■■Jn) ^ we have 

(*1 ,• • • , ■ ■ ■ ), O'l, ■ ■-jir ,• • • jn) ’ ' 

where (ii,..., i^,... in) denotes the vector (ii,... , ir-i,ir+i, ■ ■ - in) € 

This procedure defines a projection TTn of onto For each fixed r and ir € M, there 

exists a pair of matrices {Pn-i,Qn-i) which, according to formula (|3.1I1 . satisfies the equation 

= Pn-lA^^^Qn-U 

where Pn-i is a {rrp~^ x m”)-dimensional matrix with zero entries except exactly a unique entry 
1 in each row and Qn-i is an {nP x m”“^)-dimensional matrix with zero entries except, again, a 
unique 1 in each row. 

Example: For m = 2, assuming compatibility of the matrix of probability transitions in 

Af^, we calculate Pi € Af 2 x 4 and Qi G M 4 X 2 for different choices of (r,ir)- For r = 1 and P = 1 
we obtain 






f ^ 11,11 

^11,12 

^11,21 

^11,22 



/ 

1 

0 




f 1 

0 

0 

0 V 

^ 12,11 

^12,12 

^12,21 

^12,22 




0 

1 


= f ^bi 

^1,2 \ 

Vo 

1 

0 

0 J 

^ 21,11 

^21,12 

^21,21 

^21,22 




1 

0 


V ^2,1 

^2,2 J 





V ^22,11 

^22,12 

^22,21 

^22,22 

J 


V 

0 

1 

J 




For r = 1 and ir = 2: 





/ Till,11 

^11,12 

^11,21 

^11,22 



/ 

1 

0 




( 0 

0 1 

0 V 

^12,11 

^12,12 

^12,21 

^12,22 




0 

1 


( ^1,1 

^1,2 \ 

V 0 

0 0 

1 J 

^21,11 

^21,12 

^21,21 

^21,22 




1 

0 


V ^2,1 

^2,2 J 




\ ^22,11 

^22,12 

^22,21 

^22,22 

J 


V 

0 

1 

J 




Instead, for r = 2 and P = 1: 




/ Till,11 

^ 11,12 

^ 11,21 

^ 11,22 



( 

1 

0 




( 1 0 

0 

0 

^ 12,11 

^ 12,12 

^ 12,21 

^ 12,22 




1 

0 


= { ^bl 

^ 1,2 \ 

0 

0 

1 0; 

^ 21,11 

^ 21,12 

^ 21,21 

^ 21,22 




0 

1 


V ^ 2,1 

^ 2,2 / 



\ ^ 22,11 

^ 22,12 

^ 22,21 

^ 22,22 

) 


V 

0 

1 

) 




Note that in all these examples permuting simultaneously lines of Pi and columns of Qi leaves the 
product invariant. 
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For higher levels, with m > 2, fixing r = v = 1, thanks to the lexicographical order on the 
entries of the matrices we have a standard way of performing the projections of transition 
probabilities, using that: 


Pn-i — ( Idmn-i : 0 : ... : 0 


mP' 


where ‘0’ above represents the null (m” ^)-square matrices. And 


Qn—l — 


( IdjyfTl 1 ^ 


d dy^Tt — 1 


\ j 


(3.2) 


Proposition [T] implies that given a (left) eigenvector Vn £ K™'" of its projection Vn-i in 
IR™” is again an eigenvector of More precisely we have the following proposition. 

Proposition 2. Given Vn an invariant measure for a compatible Markovian chain in the product 
space M”' represented as a (row) vector in R™'", then 

V-n—l — Vn Qn—l (^■^) 

is an invariant measure in represented as a vector in R™" ^. 

Proof. Straightforward, since formula (j3.3p represents the projection (pi)* in Proposition[T} In fact, 
each column of Qn-i acts on a fixed configuration pAsi ■ ■ ■ An) S whose sum with the first 

parameter ii ranging from 1 to m gives the desired projection. □ 

As said before, matrix Qn-i in formula (13.3p is not unique and a different choice of Qn-i leads 
to a different distribution Vn-i- Nevertheless, the choice of r = 1 and A = 1 leads to the simplest 
version given by (j3.2l) . 


3.2 Back to the main example 


We go back to the example of Section [2j We apply Proposition [2] to find precisely at which level 
3 ^ n ^ 6 the bifurcation occurs given an invariant measure for the m-point motion. 

For sake of notation, we denote by Vn an invariant measure at n-points for the unperturbed 
system (e = 0), and by an invariant measure of the perturbed system (e > 0), respectively. 
We start with n = 6. We are going to compare u)) and for different initial (n = 6) invariant 
measures. First we consider the invariant measures of both systems which contain the point 123456. 
Then, in column representation we have 


/ 1123456 \ 
p _ ll24365 

1213465 

\ 1214356 /gexl 


and 


/ ll23456 \ 

1123465 
1124356 
1124365 
1213456 
1213465 
1214356 

V 1214365 / 66x1 









where the occurrence of the symbol 1 * 112 * 31415*6 ™ column vector means that the entry (^ 1 ,^ 2 , * 3 , U, 
is strictly positive while all omitted entries are zero. In addition we always assume that the distri¬ 
butions Vq and Ug are uniform (on their respective support). 

The projections of the invariant measures can easily be performed along the hrst coordinate 
(pi)* as in Proposition [21 It means that, according to formula (j3.3p . one just has to exclude the hrst 
entry of a nonzero entry (ii,..., v), 2 < r < m, in Vj, and rearrange again, if necessary, in such a 
way that the order in which they appear in the column matrix of the reduced level corresponds to 
the lexicographic order again. 

Hence we generate a sequence of vectors, which represent the invariant distributions, Vq r\ 
Ug • • • r\ and Ug rv- • • • r\ vf, where the r\ denotes the application of the procedure of the 
previous paragraph. This yields 


= Qjvl = Ql 


/ 1123456 \ 


/ 113465 \ 


/ I 3456 \ 

1124365 


114356 

^ 0 

I 3465 

1213465 


I 23456 

U 4 = 

I 4356 

\ I 214356 / 

6®xl 

\ 124365 / 

6^x1 

\ I 4365 / 


64 xl 


rv U 3 = 


/ I 356 \ 

I 365 
I 456 
V 1465 J g3xi 


^ ^^2 = ( !^^ I Ui = 

l65 ;g2xi 


/ 0 \ 
0 
0 
0 
15 
V l6 / 


For the perturbed system, following the same algorithm, we hnd the invariant measures: 


^5 ~ QI'^6 — qI 


/ ll23456 \ 


/ 113456 \ 

ll23465 


113465 

1124356 


114356 

1124365 


114365 

I213456 


I23456 

1213465 


123465 

I214356 


124356 

V 1214365 / 

6®xl 

V 124365 / 


r\Vj=Vj for all 1 < J < 3. 

This shows that the flows exhibits a 5-point bifurcation. Moreover, we shall compare invariant 
measures of the 6 point motion of and in whose supports contain the point (1, 2,1,4,1,6). 
Hence, 


/ I1214I6 \ 


/ I214I6 \ 


/ ll 416 \ 

I121315 

^ 0 

I21315 

^ 0 

I1315 

I212425 

V^ = 

112425 

rv. 114 = 

I2425 

\ I212326 / 

6 ®xl 

\ ll 2326 / 

65 x 1 

\ I2326 / 


* 5 ,*6) 
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( 1416 ^ 


/ I16 \ 


0 

1315 

0 

ll5 

0 

0 

1425 

rv U2 = 

I25 

r\vi = 

0 

V 1326 / 

63x1 

\ 126 / 

62x1 

I5 


V l6 / 

And for the perturbed system: 


4 = 


/ I121416 \ 


/ 112426 \ 


/ ll416 \ 


/ I416 

I121415 


ll2425 


ll415 


I415 

I121316 


ll2326 


I13I6 


I316 

I121315 

£ 

112325 

„ e 

I1315 

_ e 

I315 

I212426 

rvug = 

I214I6 

rv U4 = 

I2426 

rv U3 = 

I426 

I212425 


I21415 


I2425 


I425 

I212326 


I213I6 


I2326 


I326 

\ I212325 / 

66x1 

V I21315 / 

65x1 

V I2325 / 

6^x1 

V I325 


r\ 


v\ = vl 


r\ 


vl = 


This shows that the flow also exhibits a 3-point bifurcation. 


4 Degrees of freedom for fixed /c-point Markovian characteristics 

The purpose of this section is to find formulae for the dimensions of the vector space of distributions 
of i.i.d. random maps in M = {1,2,... ,m} which generate flows in M with the same prescribed 
A:-point family of transition probabilities. 

Again, as in Section 2, we use the notation for the function 

( 1 , 2 , 3 , 4 , 5 , 6 ) !-)■ 

with G M. The stochastic flow of maps (y?n)n^o in M is generated by i.i.d. random 

variables in the space of maps with the following discrete probability distribution 

m 

y= Y. (4.1) 

where 5f^^ is a Dirac measure at the mapping The non-negative coefficients G R 

are ordered lexicographically by the sub-indices. The hrst linear restriction on these coefficients 
comes from the fact that they determine the distribution of a random variable, hence 

m 

— 1- (4-2) 

We call this the 0-level restriction for the coefficients. In general, at the A:-level, for a given family 
of transition probability in A:-point motion Pui,...,uj,,vi,...,vkJ these characteristics determine linear 
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restrictions for the coefficients given by: 

where the expression {ii,... ,ik) <\ is the shorthand notation for the following vector 

, . . . , — V \, i^i+l) ■ ■ ■ iu2 — li ^2) *112+1 1 • • • 1 *11/; —1) **fc) *11/;+1) • • • j *m) • 

Obviously, the degree of freedom (dimension of subspaces which preserve the A:-point characteristics) 
is given by minus the number of linearly independent restrictions for the coefficients 
Note that, in particular, the maximal m-point transition probabilities determines uniquely the 
coefficients In fact, less than that, the information of the transition probabilities of a 

single element {ui,... ,Um) with different entries Uj’s (i.e. it does not belong to any sub-diagonal) 
is enough to determine all the m”* coefficients In other words: the information of the 

characteristics at level m means zero degree of freedom on the choice of the random mappings. 

Given the family of transition probability in A:-point motion m.-.n/;) next result determines 

the number of linearly independent restrictions for the coefficients For a hxed m € N, let 

0 < k < n < m and consider the 2-index parameter generated recursively by 

= + + 4 ) 


with Rq = 1 for all 0 ^ n ^ m. 

Theorem 2. For a finite space M = {1,2,..., m}, given a family of compatible transition proba¬ 
bilities of k-point motion with 0 ^ k ^ m, the number of linearly independent restrictions for the 
coefficients (q:ii...i^) is given by K^, defined above. In particular, Kff = m"*. 

Proof. For any 0 < A; < n < m, let denote the number of restrictions (hi. equations) at n-point 
when /c-point characteristics are given. We are going to prove by induction on k that satisfies 
the recursive equation (14.411 . 

Initially, note that when /c = 0, it means that there is no further restrictions other than equation 
()4.2p . then Rfi = 1 for all 0 < n < m. 

Assume that the formula holds for for all n € {k — l,k,..., m}. The restrictions RI^ at 

the n-level depend also on the characteristics of the (/c — 1), i.e. it is a sum of Rfi,_i plus some new 
restrictions depending exactly on characteristics at level k. This justifies the first summand on the 
right hand side of equation (j4.4jl . We just have to describe this new restrictions depending exactly 
on characteristics at level k. 

By formula (14.,411 . considering the projection at level k from level n means that there is a subset 
of positions {ui, ...,12^} C {m,..., «„} such that 


(ilv 




■■“fe 


Pui...Uk, Vl...Vk- 


(4.5) 


The number of possible subsets with k elements of {ui ,..., Un\ determines the factor (^) in equation 
(14.41) . Moreover, each parameter vi,...,Vk has m different possibilities, hence, expression (|4.5I) 
generates a block of new equations (compare with the detailed example for m = 4 in the 
Appendix). But many of this equations are linear combination of equations already counted in 
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R^_i, hence this intersection, whose dimension is given by R^_i (by definition) has to be subtracted 
in each block. This gives the recursive equation (14.4p . 

In particular, the main result of the theorem is obtained by taking n = m. Last statement is 
also trivial by taking k = n = m. □ 


Example: One easily verifies that for any m > 3 we have that > 0. It means that 

in this case one can always find different stochastic flows of mappings which differ on the 3-point 
motion, but coincide in the levels of 1-point motion and 2-point motion. For m = 3, for example 
we have Rf = 7, R 2 = 19 and R^ = 27. Hence, there exists an 8-dimensional subspace between 
3-point restrictions (of dimension R^) and 2-point restrictions (of dimension R^’^) which is given 
by distributions in the space of mappings which preserve the characteristics of the 2-point motion. 

We present a basis for this 8-dimensional space. For pairwise different i,j,k G M = {1,2,3}, 
we consider an e-perturbation of the distributions in direction of 


fijk fjjk fikk fjki T fjji T fjkk fiji T fikij 


Note that the 2-point Markovian dynamics is not affected. Varying i,j,k, we have that the 6 
possible vectors are linearly dependent (since their sum vanishes), but any choice of 5 of these 
vectors are linearly independent. The remaining three directions of the space of 2-point preserving 
dynamics can be described by the following vectors: For i,j G M, i ^ j, consider 


fiii fiij T fijj fiji T fjij fjii + fjji f: 


IJJJ- 


□ 
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Appendix A: Flows of random permutations of ,m} 

In this appendix we gather some properties and give a partial answer for the problem of the 
number of restrictions on the coefficients for flows of bijections. As before we use the same notation 
for bijective mappings fi^...im ■ ^ but here, we have that ii...im is a permutation of 

{1,2,..., m}. The stochastic flow of bijections {'Pn)n^o in M is generated by i.i.d. random variables 
in the space of permutations with the following distribution: 

(4-6) 
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The first linear restriction on these ml coefficients comes from the fact that they determine the 
distribution of a random variable, hence 


E = 1- (4-7) 

We call this the 0-level restriction for the coefficients. As before, at the A:-level, for a (compatible) 
family of transition probability in fc-point motion Pui...ui,, vi...vk, they determine linear restrictions 
for the coefficients given by; 

(il,.. 

where the sum is taken over (m — k)\ indices. As before, varying vi,... ,Vk in the expression above 
generates a block of mlj{m — k)l equations. 

In any level k, the diagonal and its complementary are invariant sets for the dynamics of random 
permutations. Moreover, for flows of bijections in a finite space, given the sub-maximal (m — 1)- 
point transition probabilities, they already determine uniquely the maximal m-point transition 
probabilities, hence they also determine the ml coefficients 

Let u = {ui,... ,Uk) and v = {vi,... ,Vk) he elements in M^. Since the order of the entries of 
the elements in does not matter in a flow, then, if cj is a permutation of k elements, then the 
transition probabilities satisfy: 


Pu^ .. .Uf^ fVl.. .Vj^ ■ ■ ■^cr(fe) VfT(l) ■ “^a(k) * 

Consider now u' = ..., and v' = (u'^, • • •, elements in such that, as 

subsets, they complement u and v respectively, i.e. {u} U {u'} = {u} U {v'} = M. Then 


E Pui-Uk,v 
ctGA 


?gA' 


(4.9) 


where A are permutations on k elements and A' are permutations in (m — k) elements. This 
is obvious from the observation that in a flow of bijections, the whole set {u} is sent to {u'} 
(independently of the order), if and only if its complementary {u} is sent to {v'}, the complementary 
of {u'}. For example: 

Pi,I = E P2...m , 5(2)^(3)...C(m)- 

5ePer({2,3,...,m}) 

The flow of bijections induced in the fc-point level sends each whole fibre (component) into a 
whole fiber. Again by the Birkhoff-von Neumann theorem the matrix of transition probabilities in 
/c-point level is again bi-stochastic for all 1 < /c < m. As we have pointed out in the introduction, in 
our context, when one deals with permutations, it means that one enters in the theory of Birkhoff 
polytopes, with many open problems. At the moment, we can give only a partial answer for the 
problem of number of restrictions on the coefficients for a flow of bijections: 

Proposition 3. For a finite space M = {1, 2,... , m}, given probability transitions of 1-point 
motion, the number of linearly independent restrictions for the coefficients (Q;q...i^) is given by 

= (m- 1)2 + 1. 
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Proof. The bi-stochastic m x m-matrix of transition probabilities of the 1-point motion has 2m — 1 
redundancies by dehnition. These redundancies corresponds to linearly dependent equations of 
type (I4.8p with k = 1. Hence the restrictions are given by [m^ — (2m — 1)] hi. equations added to 
the 0-level restriction. □ 

The arguments in the proofs of Theorem [2] and of Proposition [3] are not easily extensible to 
higher levels k in the case of bijections. This is due to the fact that, in this case, for k > 1 there 
are further restrictions which involve crossed equations coming from different blocks of equations 
generated by each fixed uiU 2 ■ ■ - Uk in equation (|4.8p (in contrast to the previous case of arbitrary 
self-maps). Moreover, for level k > m/2, new restrictions, coming from equation (14.9p which 
represents further dependence on lower levels (m — k) < k, arises (again different from the case of 
self-maps). Therefore, combinatorially, it looks non trivial to control the restrictions coming from 
different properties with non-empty intersections. As far as our knowledge, in the case of flow of 
random bijections, the problem of number of restrictions on the coefficients, given the transition 
probabilities of /c-point motion, for A: > 1, is still open. 


Appendix B: Flow of maps with m = 4 

In this appendix we illustrate in details the arguments used in the proof of Theorem [2l We consider 
here flows in M = {1,2,3,4} where there are = 256 different coefficients ai^j 2 * 3 i 4 - 

0-point motion: 


^ ^ C^ijkl — 1 
ijki 

Restrictions are given simply by Rq = 1 
1-point motions: We have (f) blocks, each block with m^ new equations: 

^ ^ C^uijk — Pl,uj U € M, 
ijk 

'y ^ (^iujk P2,ui U G Af, 
ijk 

y ^ C^ijuk — PHyUi U € Af, 
ijk 

y ^ C^ijku — PA,ui U € Af. 
ijk 

In each block we have Rq = 1 linearly dependent equations which has to be subtracted from 
the total number of equations in the block. Hence, linearly independent restrictions are given by: 
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2-point motions We have ( 2 ) blocks, each block with vn? new equations: 


E 

^uvij ~ 

~ Pl2^uv} 

u, 

V 

G 

M, 

b 







E 

^iuvj ~ 

~ P23,uv} 

u, 

V 

G 

M, 

b 







E 

^ijuv ~ 

~ P34:,UV} 

u, 

V 

G 

M, 

b 







E 

O^uivj - 

— PlS^uv-) 

u, 

V 

G 

M, 

b 







E 

^iujv ~ 

~ P2A^uvi 

u, 

V 

G 

M, 

b 







E 

^uijv ~ 

~ Pl4:^uvi 

u, 

V 

G 

M. 




In each block we have R\ = 1 linearly dependent equations (obtained by putting together the 
reduction from 2-point motion to 1-point and 0-point motion) which has to be subtracted from the 
total number of equations of the block. Hence, linearly independent restrictions are given by: 



1-M26(16 - 1 - 2 * 3) 

13-k6*9 

67. 


Remaining degrees of freedom 4^ — 67 = 189. 


3-point motions We have ( 3 ) blocks, each block with new equations: 

^ ^ O^uvwi — Pl23,uvw j UjVjlV (z M, 
i 

'y ^ O^uviw — Pl24,uvun U^V^W G iH, 
i 

y ^ Otuivw — Pl34,ni;iO) U^V^W G iH, 

j 

y ^ O^iuvw — P234,uv'w^ U^V^W G M. 


In each block we have = 37 linearly dependent equations (obtained by putting together the 
reduction from 3-point motion to 2-point, 1-point and 0-point motion) which has to be subtracted 
from the total number of equations of the block. Hence, linearly independent restrictions are given 
by: 
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= 67 + 4[64 - [10 + 3 * 9] = 175. 


4-point motions We have a single (^) = 1 block, with new equations: 

O^uvwx — Pl234:,uvwx j i U, V, W, X (z M, 


In this single block, we have = 175 linearly dependent equations (obtained by putting together 
the reduction from 4-point motion to 3-point, 2-point, 1-point and 0-point motion) which has to be 
subtracted from the total number of equations of the block. Hence, linearly independent restrictions 
are given by the following equation, which one easily sees that has a telescopic cancellation: 



= 4 \ 


We finish this appendix with the numbers for m = 5, just for comparison, without details. We 
have = 21, R^ = 181, = 821, Rl = 2101 and = 5^ = 3125. 
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